Quantum Discord: A Dynamic Approach in Geometric Picture 
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We present a dynamic approach to study the quantum discord and classical correlation. By 
local filtering operation, the evaluation of quantum discord is closely related to quantum channel 
and channel capacity. As a consequence, the traditional optimization over horizontal-or-vertical 
von Neumann measurements is replaced by that over horizontal-and-vertical three-element POVM 
measurement, from which more rigorous results of quantum discord are obtained. 
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Among various measures of quantum correlation, 
quantum discord (QD) [1] has attracted considerable at- 
tention. Many works have been devoted to quantifying 
and interpreting this measure [2-13]. However, the def- 
inition of QD involves an optimization procedure over 
the set of POVM measurements, and thus attacking the 
general case is a formidable task. Even for the general 
2-qubit states, the exact expression of QD has not been 
obtained. Restricted in the case of 2-qubit X states and 
von Neumann measurement, Ali presented an analytic 
method [10]: the direction of the optimal measurement 
is either vertical (i.e., measuring the Pauli operator az) 
or horizontal. If the X state takes the canonical form, 
the horizontal measurement is in fact measuring ax ■ Re- 
cently, it has been pointed out that there exist some ir- 
regular states for which the direction of the optimal von 
Neumann measurement is neither horizontal nor verti- 
cal, meaning that this problem has not been solved com- 
pletely even in the restrict case [14? ? ]. Meanwhile, 
there are some works devoted to finding more rigorous 
results of QD by considering more general POVM mea- 
surement [15]. 

In studying QD, the key point is to optimize the Holevo 
quantity of the postmeasurement ensemble for one sub- 
system, say system A, over all POVM measurements per- 
formed on the other subsystem B, which gives rise to the 
classical correlation. It implies that there is a close rela- 
tion between QD, in particular classical correlation and 
the channel capacity. Roughly speaking, if a channel ca- 
pacity is achieved with an ensemble containing three or 
more states [16], we can immediately conclude that there 
must be some states for which the POVM measurement 
with three or more elements is optimal to acquire the QD 
or classical correlation. In the context of 2-qubit state, 
this situation implies that the von Neumann measure- 
ment is not optimal. 

We present in this paper a dynamic method to inves- 
tigate the 2-qubit state QD. Together with the geomet- 
ric picture presented in [17], we establish an infinitely- 
many to one correspondence between the states and the 
quantum channels. The characteristic of the channel in 
turn determines whether or not the optimal measure- 



ment will allow for three or four POVM elements. In this 
framework, the traditional optimization over horizontal- 
or-vertical measurements is replaced by that over the 
horizontal-and-vertical measurements. It turns out that 
the irregular states can be explained consistently. 

Furthermore, we think that the extension of the sud- 
den transition phenomenon presented in [18] to general X 
states should be reconsidered carefully. We have found 
that for many X states the transition is smooth rather 
than sudden. Then we conjecture that, except for a spe- 
cial class of states, there is no sudden transition of QD 
for general X states undergoing two-sided phase damping 
channel. 

We consider a bipartite quantum state p^^ , and 
POVM measurement on B with a set of rank one opera- 
tion elements Mk ^ 0, J2k ^k — 1- Then QD is defined 
as Q ~ T — C. Here T is quantum mutual information, 
I = Sip"") + S{pB) - S{p^^), with p^(s) = TTBiA) P^^ 
being the reduced states and S denoting von Neumann 
entropy. The classical correlation C is given by 
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where pk = Tr[p^^(l <g) Mk)] and p^ = Tvb[p^"{1 ® 

Mk)]IPk- 

The main idea in the method is fixing the steering el- 
lipsoid and moving the reduced state. 

(i) Given a 2-qubit state p"^^, we express it as the 
unnormalized form (A ® F)p , where A is a quantum 
channel acted on qubit A, and F a filtering operation on 
B. 

(ii) The channel A determines the quantum steering 
ellipsoid <E. We prove that the form of 2; is invariant 
under any local filtering operation on B. 

(iii) Return to X state p"^^ in the canonical form. 
Having obtained the channel A associated with p^^ , we 
transform Beh state $ = |$)($| with [*) = i(|00) + |ll)) 
as $ — > [A (g) E!(^)]<I>, with S(<^) a parameterized local fil- 
tering operation. 

(iv) Denote by a"^^ {£_) the normalized form of the state 
obtained in step (iii). By selecting suitable filtering op- 
erator S(^), the a-^^{S^) remains in the X form, and 



furthermore we can translate the position of cr^{0 = 
TibO'^^{^), in the geometric picture, along vertical di- 
rection between the lower and the upper apex of the el- 
lipsoid £. 

(v) For each a^{£,), we calculate two Holevo quantities 
X^iO ^^d Xj(C)) which correspond to measuring ax and 
az, respectively, of qubit B. Compare the two Holevo 
quantities. We conclude that if Xf^iO intersects x$(0 
at a specific £,, then there are in the neighborhood of £, 
the states for which the optimal measurement is three- 
or four-element of POVM measurement. 

Now we begin with step (i). Usually we express 
the 2-qubit state p"^^ in the Hilbert-Schmidt space as 
R = 2T{p"^^)^T^ [19], where the superscript R denotes 
the reshuffling transformation: \ij){i'j'\ — > \i'i-'){j'j\, and 
the superscript T means transposition. For the sake of 
clarity, we consider X states in the canonical form. The 
canonical X state p^^ and the unitary matrix T are given 
by 
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where a + b + c + d — 1, u,v ^ 0, u^ ^ ad and v" 

Let filtering operation take the form F = 
diag(-\/a + c, \/b + d). The corresponding superoperator 
is Ap = F ® F* ^ F 1^ F. To determine the chan- 
nel A, we apply local operation A ^ Ap on Bell state 
$ and write {p^^)^ = 2A<i>-'^A|^, where the factor 2 
comes from normalization. By noting $^ — 1/2, we 
have A = {p^^)^Ap^ . Thus step (i) is finished. 

It is helpful to express the channel A in Heisenberg 
representation: 
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The effect of L on one qubit state is the transformation 
of Bloch sphere to a ellipsoidal surface, which is given by 
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where 
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FIG. 1: Schematic plot of the largest vertical intersection 
given by (xjixf' -\- [((2 — Zf^^ji^^ = 1. Two pairs of points, 
{E,F) and {G,H), correspond to the horizontal and vertical 
decomposition of p , respectively. 



It is just the quantum steering ellipsoid £ which we have 
introduced to discuss quantum discord. 

In Fig. 1, we plot schematically the largest verti- 
cal intersection of <B. The reduced state is represented 
by the point A. The ax measurement on qubit B will 
induce the (horizontal) decomposition of p^ as £^ = 
{pe,Pe'^Pf,Pf} with p^ = pp ^ 1/2 and S{p^) = 
S{pp). The (Jz measurement will give rise to the (ver- 
tical) ensemble for p^: £^ = {pg,Pc'^ Ph,Ph}- See [17] 
for detail. 

Now we prove the following lemma. 

Lemma 1. The form of the quantum steering ellipsoid 
is invariant under the loeal filtering operation on qubit 
B. 

The equation of the steering equation is given by 
{1 X y z) €{1 X y z)"^ — 0, where the symbol (£ is also 
used to denote the 4 x 4 to symmetric matrix given by 
€ = R^'^rjR-^ with 77 = diag(l, -1, -1, -1). With arbi- 
trary filtering operation S applied on qubit B, we have 

(£ ^ (£' = R-'^L^^riL^'^R-\ 

It is easy to see that L-B.rjL~ ~ |det(S)|^?7. Then it 
follows that the matrix £' is proportional to the matrix 
€ and that the equation of steering ellipsoid in invariant. 

Lemma 1 holds for any 2-qubit states. It means that 
given a steering ellipsoid £, there are infinitely many 
states corresponding to it. The reduced state of qubit 
A is represented by a point in the interior of £. We can 
move the point by performing the local filtering on qubit 
B. Step (ii) is then realized. 

We now proceed to step (iii). Define a parameter- 
ized filtering operator S(^) = diag(^, yl — C^) with 
^ G (0,1). Apply local operation A (g) A^ with A^ — 
E!(^) (g) S*(^) on Bell state $. The computable form is 
given by, after normalization. 
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The equivalent expression is 
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with i?$ = 2T$«T^, L given by (3) and L^ = TA^T^. 
When ^ = y/a + c, i?^ becomes the i? matrix of the 
X state (2). When ^ changing from to +1, the z- 
component of the Bloch vector of p^{^) changes from 
5TS ^° oTc' ■'^hich are the two apexes in vertical direc- 
tion. So step (iii) is realized. 

Step (iv) is straightforward. The steering ellipsoid (£ is 
given by (4). Any point on z axis between point G and 
point H, which represents the reduced state a [S,), is in 
one-to-one correspondence to ^ G (0, 1). The two Holevo 
quantities are given by 

^S[a^m-S[a^{0], (7) 

X$(0 - S[a^{i)] -PGS[a^{i)] -pHS[ajj{S,)] (8) 

Note that x^(0) = Xj(0) and X«-(l) — Xj(l)- It can 
be verified that there is at most one intersection point of 
Xo(0 scad xt(0 except the trivial ^ = 0, 1. 




FIG. 2: Schematic plot of X<-i-(-2) s-nd X-i{z). For any z £ 
(21,22), there exits a three- or four-element ensemble for 
a (z) which gives rise to the more optimal Holevo quantity. 



Now it is time to verify the conclusion presented in 
step (v) to complete our argument. We have to compare 
two functions x^r^{0 a-^id X$(C)- Fo^' the sake of clarity, 
we replace the argument of the Holevo quantities by the 
position of cr (C) in the geometric picture, denoted by z. 

We plot schematically in Fig. 2 the functions x+^{z) 
and xxi^)- Assume that they intersect at z. We plot a 
line tangent to both curves with tangent points denoted 
by M and N . Corresponding to M and N , the reduced 
states are (t^{zi) and a^{z2), respectively. It can be seen 
that z G (zi,Z2). For a (zi), the horizontal decompo- 
sition, £^{zi) = {pE,cr^izi); pf,<Jf{zi)}, gives rise to 
X<->.(-Zi)- Here, both probabilities pe and pp are in fact 
equal to 1/2. For a^{z2), the vertical decomposition, 
Sf{z2) = {PG, 0-0(^2); Pff,cr^(z2)}, gives rise to x$(^2)- 

Now consider a state a-^{zq) = qa'^{zi) + {l—q)a^{z2) 
with q G (0, 1). We can see that Zq = qzi + (1 — q)z2- Let 



us consider the 4-state ensemble for a^{zq): 

^iizq) =[<lPE,cri{zi); qpF,cr^izi); 

(1-9)^0,^0(^2); {l~q)PH,(THiz2)}- (9) 

We calculate the Holevo quantity of the ensemble £f(zq). 

X4{zq) = qx^{zi) + (1 - Q)xt{z2) + S[a^{zq)] 

-qS[a''{z,)]-{l-q)S[a^iz2)]. (10) 

The linear combination of the last three terms at the 
right-hand side of (10) is strictly larger than zero due 
to the concavity of entropy. As for the first two terms, 
it is easy to see that for any q G (zi,Z2) the sum 
QX^{zi) + {l-q)xt{z2) is strictly larger than both x^{q) 
and x^il)- It turns out that Xi{zq) > X-^il) and 
Xi{zq) > X^il) f"^^' 9 G (zi,Z2). Thus we prove that, 
for states in this interval, four-element POVM measure- 
ment is optimal. In fact, four-element ensemble £4 can 
be simplified to three-element ensemble £^^ in which one 
element corresponds to one of two the vertical apexes of 
€, and the other two elements correspond to two intersec- 
tion points of a vertical line and £. From this perspective, 
we conjecture that the Holevo quantity x{z) could be a 
smooth curve. 

In the context of evaluating QD or classical correlation 
of X states, the optimal three-element POVM measure- 
ment is such that one measurement operator is |0)(0| or 
|1)(1| and the other two are | + x){+x\ and | — x){~x\ 
respectively. Applying this type of measurement to the 
irregular state presented in [14], we see a more rigorous 
result of QD. 

The dynamic method presented above is concerned 
with the situation that the quantum steering ellipsoid 
is fixed while the reduced state of one qubit is pushed up 
and down by the filtering operation applied on the other 
qubit. On the contrary, we can fix the reduced state 
and deform the ellipsoid. This process can be realized by 
using phase damping channel. 

When two-sided phase damping channel is acted on 
the state p"^^ , the time evolution is given by 
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where Ki — diag(l,7) and K2 — diag(0, \/l — 7^) are 
Kraus operators representing phase damping channel, 
and 7 = e^'"* with F the phase damping rate. Here 
we assume that qubits A and B endure the same noisy 
environment. At initial time t = the steering ellip- 
soid (£(0) is given by (4). At time i > 0, the ellipsoid is 
transformed to €(i), which is expressed by 
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That is, with 7 decreasing from 1 to 0, the radius of 
the elhpsoid along x axis and that along y axis decrease 
continuously from li and I2 respectively to zero, whereas 
the radius along z axis remains the same. The reduced 
state p^ is not affected by phase damping. For the state 
p^^{t) at time t, we have two Holevo quantities: x^{t) 
and xt- Note that x± is independent of time t. If initially 
X-i^iO) > X$i there must be a critical time i on which 
X<->.(f) = Xj- It follows that x^{t) < Xj for all t > t. 
Horizontal-or- vertical optimization tells us that at time i 
there is a sudden transition from classical to quantum de- 
coherence regime. However, the horizontal-and-vertical 
optimization over 3-element POVM measurements shows 
that this sudden transition may not exist. 

For example, let us consider the state presented in [20]. 
We compare in Fig. 3 two numerical results of the clas- 
sical of the state undergoing two-sided phase damping 
channel with F = 0.01. It is shown that the three-element 
POVM measurement is more optimal than the traditional 
von Neumann measurement. 



As a consequence, the traditional optimization over von 
Neumann measurements is shown to be not sufficient to 
give rigorous results. 

This work was supported by National Nature Science 
Foundation of China, the CAS, and the National Funda- 
mental Research Program 2007CB925200. 
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FIG. 3: Time evolution of the classical correlation of the 2- 
qubit state, with the element of p ^ given by a = 0.4875, 
b = 0.1625, c = 0.0875, d = 0.2625, u = 0.3354 and 
V = 0.1118, undergoing two-sided phase damping channel. 
Solid line corresponds to the result coming from horizontal- 
or- vertical von Neumann measurement, while the dotted line 
to the result obtained by three-element POVM measurement. 



A special case must be pointed out. If the center of the 
ellipsoid £ is on the origin point, the sudden transition 
must occur as long as x<-i-(0) > X$- In this case, when 
X<->.(i) is equal to xj at the critical time i, the largest 
vertical intersection of the €(i) is a circle. It turns out 
that X<->(^) ^) is identical x$(-z) for all z between the lower 
and the upper apexes of the circle, and there is no region 
in which the four- or three-state ensemble can yields a 
more optimal Holevo quantity. 

In conclusion, we present a dynamic approach to study 
quantum correlation and classical correlation. With the 
the dynamic process illustrated explicitly in the geomet- 
ric picture, we find that the QD can be studied from 
the viewpoint of quantum channel and channel capacity. 
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